Heavy-light meson decay constants with Nf = 3 
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During the past year the MILC Collaboration has continued its study of heavy-light meson decay constants 
with three dynamical quarks. Calculations have been extended to a second lattice spacing of about 0.09 fm. At 
this lattice spacing, there are results in the quenched approximation and for three sets of dynamical quark mass: 
mi = m s ; mi = 0Am s and m; = 0.2m s , where mi is the light mass for the u and d quarks and m s is the strange 
quark mass. At the coarser lattice spacing, for which results were presented at Lattice 2001, statistics have been 
increased for two sets of quark masses and three additional sets of quark masses have been studied, giving a total 
of eight combinations used to interpolate between the quenched and chiral limits. When these calculations are 
completed, we can study the decay constants taking into account both chiral and continuum extrapolations. 



1. INTRODUCTION 

We are extending a calculation of heavy-light 
meson decay constants with three flavors of dy- 
namical quarks that was begun last year fjl]. In 
addition to increasing statistics on some runs, we 
have additional mass combinations for a = 0.13 
fm and new dynamical quark runs with a = 0.09 
fm. Runs for the coarser lattice spacing are com- 
pleted, but running will continue for the finer lat- 
tice spacing, so that we can begin to understand 
the continuum limit. 

Dynamical gauge configurations are generated 
using the Asqtad action pi. For the heavy- light 
mesons, we use tadpole-improved clover valence 
quarks and operators that are improved accord- 
ing to the Fermilab formalism For each en- 
semble of dynamical quark configurations, we use 
five light and five heavy valence quark masses. 

* presented by S. Gottlieb 



The masses and decay constants are interpolated 
or extrapolated as explained below to get physi- 
cally relevant values with the overall scale set by 
the p mass. 

2. PROGRESS SINCE LATTICE 01 

Last year we had completed running for a — 
0.13 fm on quenched configurations, on a two fla- 
vor ensemble, and for four combinations of light 
and strange quark dynamical masses. A fraction 
of another set was run. We have now completed 
eight sets of dynamical masses with a — 0.13 fm. 
Our quenched run for a = 0.09 fm included about 
15% of the configurations there, and we were gen- 
erating configurations with two sets of dynamical 
masses. We soon completed analysis of every- 
other quenched lattice and have made substan- 
tial progress on three dynamical ensembles. The 
table below contains a summary of our running. 
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Figure 1. Crural extrapolation of fs showing 

both linear and non-linear fits for (3 — 6.79, 
kq = 0.08. 



3. ANALYSIS OF RESULTS 

The analysis of the heavy-light decay constants 
involves a number of steps: 

1. fit light pseudoscalar masses 

2. perform chiral fit of pseudoscalar masses to 
determine k c 

3. fit light vector meson masses 

4. perform chiral fit of vector meson masses and 
determine m u ^ and a to get physical jm p and 
m p 

5. determine m s from mass of ss pseudoscalar 
state assuming linear chiral mass relation 

6. fit heavy-light channels to determine masses 
and decay amplitudes 

7. extrapolate or interpolate results in light 
quark mass to m uc i or m s , respectively (see 
Fig. 1). One can see that the confidence level 
of the nonlinear fit is much better than that of 
the linear fit. The curvature of the former has 
the same sign as that expected from chiral log- 
arithms as recently suggested 0; however, our 
data has considerably less curvature than the low- 
est order chiral logarithms would predict, requir- 
ing higher-order terms and, perhaps, a modifica- 



tion of the logarithms due to 'taste' breaking or 
partial quenching. 

8. after removal of perturbative logarithms, fit 
fQq\Z m Qq to a power series in l/rriQ q and inter- 
polate to B, B s , D and D s meson masses 

9. put the perturbative logarithm back and 
use the heavy-light axial- vector current renormal- 
ization constant to get the renormalized decay- 
constant 

Unfortunately, the axial- vector renormalization 
constant has not yet been calculated either per- 
turbatively or nonpcrturbatively. We assume 
that in the quenched approximation our current 
results at a — 0.13 fm with an improved action 
agree with the continuum limit of our previous 
calculation using Wilson and clover quarks and 
the Wilson gauge action. This was explained in 
more detail in Ref. ||. 

After the above procedure is done on each en- 
semble, we have a partially quenched result at a 
particular value of dynamical m^/mp. We then 
plot these results as a function of (m v /m p ) 2 to 
perform a chiral extrapolation. This is demon- 
strated for Jb in Fig. 2. 

It is worthwhile to plot the ratio of B s and B 
meson decay constants since many of the system- 
atic errors are common, and a good deal of the 
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Figure 2. f B as a function of (m^/mp) 2 
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4. CONCLUSIONS 
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Our preliminary three flavor results and er- 
ror estimates are slightly changed from last year. 
With higher statistics the quenched results at the 
finer lattice spacing are in good agreement with 
the coarser spacing for the B mesons. We find: 



W/B UenCh [ m p scale ] 



1.23(4)(6), 



/ s .// fl = 1.18(1X1?). 

If we take the MILC continuum quenched value 
(m p scale) for Jb (169 MeV) as given, our value 
/s//s UCnch = 1-23 gives f B « 207 MeV with 
2+1 flavors. An independent determination of fg 
awaits the calculation of the axial current renor- 
malization factor. The errors are detailed in the 
table below. 

Better understanding of the chiral logs is neces- 
sary and may change the fg chiral extrapolation 
beyond the error estimated here. 
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